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We describe the collapse of the bosonic component in a boson-fermion mixture due to the pressure 
exerted on them by a large fermionic component, leading to collapse in a system with all-repulsive 
interactions. We describe the phenomena early collapse and of super-slow collapse of the mixture. 
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Introduction.- The achievement of Bose-Einstein con- 
densates (BECs) Q has pushed the field of degenerate 
quantum gases to become one of the most active areas 
of Physics. After the condensation of different bosonic 
species, degenerate mixtures of bosonic and fermionic 
atoms were created providing a highly controllable 
tool for study of systems of mixed quantum statistics. 

Several interesting nonlinear phenomena have been 
studied in the context of boson-fermion mixtures. First, 
the fact that interspecies interactions may result in at- 
traction among bosons |3|] , implies that solitons could also 
in quasi-one dimensional boson-fermion mixtures 
Also, collapse of the atomic cloud induced by the in- 
terspecies attraction in boson-fermion mixtures was ob- 
served experimentally 0] and studied theoretically 

The usual scenario of collapse in single-species quan- 
tum gases corresponds to a BEG which collapses because 
of its attractive self-interaction 0. This phenomenon 
had been long known in mathematical physics but 
the theoretical description of its specificities for BECs has 
motivated a lot of theoretical research. In boson-fermion 
mixtures, it is the attractive interspecies interaction what 
drives the blow-up of the atomic cloud 

In this paper we describe new blow-up scenarios in de- 
generate quantum gases, namely in boson-fermion mix- 
tures, when all interatomic interactions are repulsive. 
More specifically we will describe the phenomena of early 
collapse and of super-slow collapse of the mixtures. 

The model.- We describe bosons in the mean field ap- 
proximation and spin polarized fermions at zero tempera- 
ture in the hydrodynamic approximation, when the num- 
ber of fermions is much larger than the number of bosons. 
The respective dynamical system reads 0, [Tli 



2mh 



A^-f H(r)V' + g&f,|'0lV + 5&/P'0, (la) 



V 



. 1/3 + — IV'I^ 



(lb) 



where 'ip{r,t) is the macroscopic wave function describ- 
ing the bosonic component and pi{r,t) is the excitation 



of the density of fermions, due to boson- fermions interac- 
tions, from the unperturbed density poir), so that the to- 
tal fermionic density is given by p{r, t) = po{r)-\- pi{r, t). 
gtb and gb / are the coefficients of two-body boson-boson 
and boson-fermion interactions, Vb(r) is the confining po- 
tential for bosons, V& = mba;^r^/2. 

It is convenient to introduce the scaled variables 
X — r/a and r = ajf,^, where a = ^ h/mbWb is the 
linear oscillator length for the bosons. We also de- 
fi ne Tip ^ {gbf /fihJb)po, n = [gbf/hhJb) pi and u = 

gbfmb/huJbmfTp, and rewrite Eqs. ([T]) in the form 



--V'^u + Vou + nu + g\u\'^u, (2a) 



UT-r = V 



(2b) 



where Vq = Vo{x) = -\- no{x), g = gbbrnf/gbfirib, 

a = (67r^)^^^ h^/^mb/'im?^ (gbf)"^^^ ^l^^y ^ is understood 
in terms of the new coordinates and x = |a;| for a; € K'^. 
We will refer to this model as the nonlinear Zakharov 
system (NZS) since it can be seen as a generalization of 
the Zakharov system [l^ describing Langmuir waves in 
plasmas, the latter being given by Eq. ((2]) with gbb — 
(i.e. in the absence of the the two-body interactions 
among bososns). We will consider Eqs. ^ with with 
analytical initial data u{0,x), n{Q,x), and nT-{0,x). 

Conserved quantities.- Equations ([2]) can be rewritten 
in a Hamiltonian form. To this end, following [l3|, we 
introduce the vector function 



(3) 



V = - / TioV ( —n + \u\^ ) dr, 



which describes the hydrodynamic velocity of excitations 
and rewrite the NZS in the form 



= — Vv, = — tiqV 



1/3 



(4a) 
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which admits the Hamiltonian 



H 



I' + Vilup 



1 n 



1/3 



2 a 



^n+\u\ 



da;. (5) 



(hereafter all integrals without explicit limits are taken 
over the whole space). H is an integral of motion: 

Another integral of motion is given by the total num- 
ber of bosons, Nb = \\u\\2, defined through the standard 

notation for the norm: ||u||p = (/ \u\PdxY^^ where p 
is an integer. We also notice the n describes excitations 
of the Fermi sea and thus the number of excited fermions 
is not conserved. 

Modulational instability.- The simplest manifestation 
of instabilities in nonlinear wave systems is the appear- 
ance of modulational instability. Neglecting the external 
trapping potential, i.e. supposing that no, and hence 
go = nQ^^/a, are ^-independent constants, and look- 
ing for plane-wave solutions of Eq. ([2a|) of the form 
Up = ^exp(iKa; — iflr) we get the corresponding ex- 
citation of the fermionic density Up = —goA^. Next we 
look for perturbed solutions of the form u — Up{l + ui), 
n = np(l + rti), where ^ \up\ and |ni| <^ |rtp|. After 
linearizing with respect to (ui,ni) oc exp(zka; — a;i), from 
the resulting linear system we find the dispersion relation 

2A^gok'^ = (cj2 _ j^2^^^2 _ . j^)2 _ 2gA^e - k^] (6) 

The instability occurs when at least one of the frequencies 
Lo obtained from Eq. ^ has a nonzero imaginary part. 
This happens for the wavevectors satisfying 



< 2A^ {go -g)- 4:K^ cos^ 0, 



(7) 



where 6 is the angle between K and k. Thus the condi- 
tion of modulational instability has the same functional 
form as that for the purely bosonic condensate [ifli] with 
the substitution of —g by go — g. i.e. if go > g, then 
even when all two-body interactions are repulsive pro- 
vided g < go modulational instability may set in. From 
the physical point of view, the modulational instability 
develops for long wavelength excitations. For such excita- 
tions the fermionic excitations are smooth, and therefore 
the constant amplitude density Up enters the condition 
of the instability. Thus, this result means that instabili- 
ties may develop in a system with repulsive interatomic 



interactions, which is very different from the case of the 
purely bosonic condensate and is the first indication of 
the phenomenona to be discussed in detail in this paper. 

Non-collapsing solutions.- First, we will prove that 
there are initial data which do not blow, i.e. the existence 
of global solutions of Eq. More specifically, we want 
to establish imiform boundness in time oi U = ||Vm||2 
13, lJ| for g < go (for g > go the proof is trivial). Using 
the simplified version of the Gagliardo-Nirenberg inequal- 
ity ||u>||| < Ci||Vw||^|lw||^ ilj, with Ci = 0.44927... 
being the best constant in the 3D case [l5|, we get 



(8) 



1 /2 

where 7 = (50 - g)CiN^' . Eq. ([8]) is vahd for aU r, 
and thus, if there is at least one positive root of F{U) = 
-fU^ - -f = 0, to be denoted as u+: F{u+) = 0, 
and if initially U{Q) < u+, then Eq. ([8]) guarantees that 
U (t) < M-i- for any time, i.e. U is bounded uniformly in 
time. Hence, to establish the sufficient condition for the 
existence of the solution we have to find the condition 
of the existence of a positive root u+. It is easy to see 
that F{U) has a local minimum at Umm = 2/87 and a 
sufficient condition for the existence of a positive root is 
F{Umin) < 0, what is equivalent to the constraint 



H<Ho = 



4 1 



2772 27 (go - g^CfNt, 



(9) 



It can be seen that the condition U{0) < u+ is satisfied by 
choosing U{0)^ < H. This rigorous result, being a suffi- 
cient condition, provides a lower bound for the largest en- 
ergy allowing the existence of solutions, i.e. the physical 
region of existence could be larger, in particular allowing 
H > Hq. Since Eq. © involves both the Hamiltonian 
H and the number of bosons A^;,, which are not indepen- 
dent, in Fig. [T]we plot the region of non collapsing ground 
states by fixing the number of fermions and computing 
numerically the ground state for different A'^;, until the 
limit of the inequality is reached. This leads a domain in 
which we can guarantee the existence of global solutions. 
It is interesting to note that there exists a critical number 
of fermions, A^* below which all solutions exist globally 
and no collapse can occur. The number of fermions sepa- 
rating the domain of globally existing solution decreases 
with Nh, what is explained by the fact that the repulsion 
between bosons due to two-body interactions decreases 
with Nb, and thus smaller fermionic densities are able to 
dominate the inter-bosonic repulsion. 

Collapsing solutions.- It was shown numerically in 24 1 
that the Zakharov equations may undergo finite-time col- 
lapse. In our system, due to the presence of repulsive 
boson-boson interactions it may seem counter-intuitive 
that such a phenomenon could exist. However, we will 
argue that the bosonic component in a Bose-Fermi mix- 
ture, in which there are many more fermions than bosons, 
can undergo collapse in finite time even when both the 
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FIG. 1: Ground state solutions of Eq. ^ for atf — 10 nm, 
att = 5.25 nm, m/ = 0.66 x 10"^^ kg, mt = 1.44 x lO'^^kg, 
uib = 188 Ifz, and cuf = LUb\/mb/mf — 269 Hz, corresponding 
to ^^Kh-^'^K mixtures such as those of Ref. 01 . Shown are 
(a) the location on the plane A^;, — Nf of the ground states 
satisfying H = 4/277'^, thus all ground state solutions in the 
shaded region must be stable. Panel (b) shows some more 
data on a logarithmic scale. 



self-interaction between bosons and the interspecies in- 
teractions between bosons and fermions are repulsive. To 
support this statement we will use "early collapse" type 
arguments [l^, modified for the NZS since the standard 
arguments based on the virial identities cannot be used. 
Let us define the squared wavepacket width as Y{t) — 



dx. By direct derivation we obtain Y{t) — 



4 Im J u{x ■ Vu) dx and 

Y{t)= I (8|Vup da; + 6g|u|''-4|wpa;Vn)da;, (10) 



(hereafter an overdot stands for the derivative with re- 
spect to r) . In our case and even in the case of classical 
Zakharov system [3| the lack of control on the last term 
of Eq. (fTO|) makes impossible to apply the usual virial- 
type arguments. Nevertheless by analyzing Eq. (fTO|) for 
small times we can extract information about the early 
stage of the evolution. First, we notice that Y(Q) and 
y(0) depend only on the initial distribution of the boson 
component, u{Q,x). In particular, for simplification of 
the further analysis, we can make 1^(0) = by choosing 
m(0, x) real. Considering the higher order derivatives one 
can note that for given u(0,£c), F(0) can be made neg- 
ative by an appropriate choice of n(0, a;). For example, 
this occurs if n(0, x) is a cup-like profile. Moreover, from 
the expression for the third derivative 



y(0) = -4 / \u{Q,x)\^x -Vuri^.x) dx, 



(11) 



it follows that y'(O) can also be made negative by an 
appropriate choice of nT-(0, x). 

Summarizing the above arguments, assuming that 
Y{t) is sufficiently smooth and requiring F(0) and y(0) 
to be negative, what is always possible by means of the 
proper choice of the initial conditions, we can expand 



Y{t) in Taylor series for sufficiently small r: 

Y{r) ^ (r^ - r,2)M + yy(0) + o{r% (12) 



Here r* = ^ -2y(0)/F(0). Now we notice that the 
higher terms that are not written explicitly in (|12p de- 
pend on higher derivatives of the fields it(r, a;) and n(r, a;) 
at r = [this can be verified by direct computation of 
the higher derivative of Y{t), by analogy with (fTU]) and 
pip ], and thus can be chosen to satisfy a priori given 
properties (which in our case are reduced to existence 
of a nonzero radius of convergence of the Taylor series). 
This means that one can construct initial data such that 
Y{t) crosses zero at r < r*. This contradiction with the 
definition of Y(t), which is nonnegative, implies a finite 
time collapse for the chosen initial data. 

The above arguments, would become mathematically 
rigorous only if one could prove that r* is smaller than 
the convergence radius of the series (IT^ . Although that 
proof is not available, we can show that by proper choice 
of the initial conditions one can make as small as nec- 
essary, what supports our conjecture. 

Let us consider a realistic situation, where initially the 
fermions are unperturbed and the bosons have a Gaus- 
sian distribution narrower than the size of the fermionic 
cloud. In that situation, no can be taken as constant and 
we can compute r*, given below in the physical units: 



(2^) 



l/4„l/2 



mj ( {QpoY/^abf 2abb 



oJbiNbObfy^'^ nib + ruf 



.1/3 



-1/2 



(13) 

where uJb is the initial width of the bosonic cloud. From 
Eq. it can be clearly seen that controlling Ob / allows 
one to make T* as small as necessary. Taking the same 
data as in Fig.[l] but with Obf — 10~^to, we get a blow-up 
time r* w 0.127 s. 

Super-slow collapse.- Finally, following Ref. [l^ where 
the possibility of collapse was studied for the conven- 
tional Zakharov system, we outline the proof of the fact 
that the bosonic component for g < go and H < may 
undergo a completely different physical phenomenon: in- 
finite time blow-up, which we will denote as super slow 
blow-up. We restrict our analysis to radially symmetric 
solutions and for the sake of simplicity consider a spa- 
tially homogeneous mixture with a = 1 and no = 1 (i.e. 
go — 1), what can always be done by proper rescaling of 
the independent variables. Recalling that x is a dimen- 
sionless radial variable we define 

poo 

Z{t) = Im / {x^uUx + xnw ■ x) fxdx, (14) 
Jo 

where the weight function f{x) is chosen, such that for 
arbitrary positive constants Ro, C2 and C3, the following 
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conditions are satisfied [23|: 

\M<Ro, <l,and |A2/|<C2, 

C3 > 3 - A^/ > for a; > i?o (15) 

3 - Arf = for a; < i?o, 

where = + § ^ • direct algebra we obtain 



-Z(r) 



\vr 2 
^A2/+(3-2/,,)K|2 + -/,|v|2 
2 X 



+ (3-A./)(nH^ + |H4 + ^^^^ 



and unisng properties (|15p we estimate 



X da; 



377 

4^' 



((3-A,/)lu|4 + |u,P)a;2dx 



ft-^. (16) 



Radially symmetric functions u G L^fR^), with S 
L^(R^), satisfy the Strauss lemma [23|: jwp < 
/3i?~^|jM||2||w2;||2 for X > R, where R is an arbitrary 
positive constant and the constant /3 is independent on 
R. Since g < 1, application of the Strauss lemma with 
R = Rq together with (fT5|) enables us to estimate the 
right hand side of Eq. (fT6| as 



Z > 



3H 

An 



^1 
2 ' 



16 R^ 



(17) 



Thus if H is negative, then choosing the constants C2 
and Rq one can obtain ~Z{t) > 0. On the other 
hand, following [l^, it can be shown that —Z{t) < 



^dx] where C4 is a posi- 



tive constant depending only on ||u||2 and on \ fx 
leads to 



lim 



^dx = 



This 



(18) 



which implies collapse for infinite times and physically 
gives an indication that a large fermionic component is 
able to slowly compress a bosonic component even when 
two-body interactions of bosons are repulsing, provided 
they are not too strong. This leads to a super-slow col- 
lapse. 

Conclusions.- We have studied different collapse sce- 
narios of the bosonic component in boson-fermion mix- 
ture at limited positive boson-boson scattering lengths 
(defined by 17 < go) provided the fermionic component 
has the density large enough. They are (i) a finite-time 
collapse, occurring at early stages of the evolution, and 
(ii) super-slow collapse happening at infinite times. The 
physical reason for collapse in both cases is the pres- 
sure exerted on the bosons by the fermionic component. 
It is remarkable that these phenomena occur even when 



all interatomic interactions are repulsive. Our theoreti- 
cal predictions can be tested with current experimental 
setups and describe novel interesting singular quantum 
phenomena with degenerate quantum gases. 
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